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ABSTRACT 


The spatial and temporal invariance in the spectra of energetic particles in the grad- 
ual solar events is reproduced in the simulations. Based on a numerical solution of the 
focused transport equation, we obtain the intensity time profiles of solar energetic 
particles (SEPs) accelerated by an interplanetary shock in the three-dimensional inter- 
planetary space. The shock is treated as a moving source of energetic particles with a 
distribution function. The time profiles of particle flux with different energies are cal- 
culated in the ecliptic at 1 AU. According to our model, we find that shock acceleration 
strength, parallel diffusion and adiabatic cooling are the main factors in forming the 
spatial invariance in SEP spectra, and perpendicular diffusion is a secondary factor. In 
addition, the temporal invariance in SEP spectra is mainly due to the effect of adiabatic 
cooling. Furthermore, a spectra invariant region, which agrees with observations but 
is different than the one suggested by Reames and co-workers, is proposed based on 


our simulations. 


Subject headings: Sun: activity — Sun: coronal mass ejections (CMEs) — Sun: particle 


emission 


zug 
1. INTRODUCTION 


Solar energetic particle (SEP) events can roughly be divided into two categories: impulsive 
events and gradual events (Reames 1995, 1999). The impulsive events, with the characteristics 
of low intensity and short duration, are produced by solar flares. Gradual events, usually lasting 
longer and having high intensity, are related to the shocks driven by interplanetary coronal mass 
ejections (ICMEs). Lario et al. (2006) investigated the radial and longitudinal dependence of 4— 13 
and 27 — 37 MeV proton peak intensities and fluences measured within 1 AU. They found the 
peak intensities and fluences of SEP events can be approximated by j = jor * exp |-«(6 as $o], 
where j is either the peak intensity or the fluence, r is the radial distance of the spacecraft, ¢ is the 
longitudinal angular distance between the footpoint of the observer’s field line and the region of 
the SEP source, and 4» is the centroid of the distributions. Furthermore, the radial dependence of 
peak intensities and fluences of SEP events have been simulated with a focused-diffusion transport 


equation (Lario et al. 2007). 


Generally, there are two major approaches to modeling SEP acceleration by CME driven 
shocks: some authors (Heras et al. 1992, 1995; Kallenrode & Wibberenz 1997; Lario et al. 1998; 
Kallenrode 2001; Ng et al. 2003; Wang et al. 2012; Qin et al. 2013) adopted a “black box” 
model to treat the shock as a moving source, and SEPs are injected at the shock with an assumed 
injection strength, while a few other studies include the acceleration of SEPs by shocks (Lee 
1983; Gordon et al. 1999; Zank et al. 2000; Li et al. 2003; Rice et al. 2003; Sokolov et al. 2004; 
Li et al. 2005; Kóta et al. 2005; Tylka & Lee 2006; Zuo et al. 2011; Zuo et al. 2013). In these 
models, three important effects of acceleration and propagation mechanisms have been involved. 
The first effect is the acceleration process by the CME-driven shock. Zank et al. (2000) modeled 
the evolution of a CME-driven shock based on an “onion shell" model, and this model has been 
furtherly developed in a number of papers (Li et al. 2003; Rice et al. 2003; Li et al. 2005). They 


used a magnetohydrodynamics (MHD) code to describe the evolution of the CME-driven shock 
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in the interplanetary space, and wave excitation by streaming energetic particles produced at 
shock is included. Based on the model, the simulation results can successfully explain the SEP 
fluxes and spectra in some multi-spacecraft observed events (Verkhoglyadova et al. 2009, 2010). 
The second effect is energetic particles to interact with Alfven wave self-consistently. Ng et al. 
(2003, 2012) presented a model of particle transport including streaming proton-generated Alfven 
waves, and the amplification of the Alfven waves is determined by the anisotropy of particles. The 
particle diffusion coefficients can be calculated from wave intensity and wave growth rates. Their 
simulation results show a good agreement with the observed spectral slope and abundance ratios 
of heavy ions. The third effect is the realistic geometry of CME and its shock (Sokolov et al. 2004; 
Kóta et al. 2005). Sokolov et al. (2004) modeled particle acceleration and transport as CME driven 
shock wave propagating from 4 to 30 solar radius from the Sun. The realistic structures of CME 
and its shock are derived from a numerical solution of a fully three-dimensional MHD model. 
Their simulation results demonstrate that the diffusive shock acceleration theory can account for 


the increase of hundreds of MeV protons during the early stages of CME driven shock. 


Solar energetic particle events measured by multi-spacecraft help us to understand the 
processes of particle acceleration and transport in the heliosphere. In some gradual events, the 
SEP fluxes measured by widely separated spacecraft present similar intensities within a small 
~ 2 — 3 factor in different latitudes, longitudes or radius (Reames et al. 1997; Reames 2010, 
2013; McKibben et al. 2001b; Maclennan et al. 2001; Lario et al. 2003; Tan et al. 2009). This 
phenomenon was firstly proposed by McKibben (1972), and was named "reservoir" by Roelof 
et al. (1992). In order to interpret the reservoir phenomenon, McKibben (1972) and McKibben 
et al. (2001b) involved an effective perpendicular diffusion to reduce the spatial gradients of 
flux, while Roelof et al. (1992) suggested a diffusion barrier produced by ICMES or shocks. The 
magnitude of magnetic field increases at the outer boundary of reservoirs, so that SEPs could 
be contained in the reservoirs for a long time. Furthermore, the interplanetary magnetic field 


(IMF) has been disturbed by ICMEs, SEPs could be redistributed. Reames et al. (1996) shows 
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that, in some gradual SEP events, the spectra are invariant both in space and time. This discovery 
extended the original work of McKibben (1972). In Reames et al. (1996), they considered an 
expanding magnetic bottle of quasi-trapped particles between an ICME driven shock and the 
Sun. As the magnetic bottle expanding, the SEP fluxes gradually decrease as a result of parallel 
diffusion and adiabatic cooling. In this sense, the magnetic bottle plays a pivotal role in the decay 


phase of SEP event. 


In principle, the disturbances in the magnetic field caused by ICMEs can help the particles 
redistribute in space. However, the reservoir phenomenon cannot be simply explained as a result 
of the disturbances of IMF caused by ICMEs. Firstly, in the redistribution process in Reames et al. 
(1996) and Reames (2013), no explicit transport mechanism can reduce latitudinal, longitudinal, 
and radial gradients of SEP fluxes besides perpendicular diffusion. Secondly, in some SEP events, 
ICMEs are not directly observed by the spacecraft, the reservoir phenomenon is also observed 
(McKibben et al. 2001a). Thirdly, in Reames (1999), when the observer is located at the eastern 
of the shock, the onset time of temporal invariance in the SEP spectra is earlier than the shock 


arrival time. These results are not consistent with that of an expanding magnetic bottle. 


The effect of perpendicular diffusion is important in the SEP fluxes especially when the 
observer is disconnected from the shock by IMF. During the time period March 1, 1979 to March 
11, 1979, a gradual SEP event has been detected by Helios 1, Helios 2, and IMP 8. The three 
spacecraft are located in the ecliptic near 1 AU, but at different longitudes. In the decay phase of 
this SEP event, the reservoir phenomenon appeared (Reames et al. 1997; Reames 1999, 2010, 
2013). In this event, the in-situ observation shows that an ICME was detected by Helios 1, but not 
by Helios 2 and IMP 8. And the interplanetary shock was only observed by Helios 1 and 2, but 
not by /MP 8 (Lario et al. 2006; Reames 2010). According to the location of the three spacecraft, 
if the ICME was located behind the center of the shock front, then Helios 1 was located near to the 


center of shock, and Helios 2 and IMP 8 were located at the West frank of the shock. However, 


cc 


the onset time of SEP fluxes observed by three spacecraft was very close. How could SEPs be 
detected by JM P 8 before it was connected to the shock by IMF? One possible answer is the effect 


of perpendicular diffusion which also possibly works in forming the reservoir phenomenon. 


However, perpendicular diffusion has been always a difficult problem for several decades. 
Observation results show various levels of perpendicular diffusion coefficients for different SEP 
events. For example, ‘dropout’ phenomenon in the impulsive SEP event (Mazur et al. 2000) 
usually show reduced perpendicular diffusion, in order to reproduce the ‘dropout’ phenomenon 
in simulations, perpendicular diffusion coefficient x, should be several order magnitude smaller 
than parallel one xı (Giacalone et al. 2000; Guo & Giacalone 2014; Dröge et al. 2010; Wang 
et al. 2014). On the other hand, for some events, observation results show that the perpendicular 
diffusion coefficients could be comparable to the parallel ones (Dwyer et al. 1997; Zhang et al. 
2003; Dresing et al. 2012). In order to understand diffusion, many efforts were made theoretically. 
By assuming energetic particles’ perpendicular and parallel diffusion do not interaction, Jokipii 
(1966) developed the quasi-linear theory (QLT). According to QLT, perpendicular diffusion 
coefficient is usually much smaller than the parallel one. However, it is found that interaction 
between parallel and perpendicular diffusion is important in theory (Kóta & Jokipii 2000) and in 
simulations (Qin et al. 2002b,a), so the non-linear guiding center (NLGC) theory (Matthaeus et al. 
2003) is developed to describe perpendicular diffusion with the influence of parallel diffusion, 
which agrees with simulations much better than QLT. In addition, simulations show different 
levels of perpendicular diffusion, e.g., in Qin & Shalchi (2012), the magnitude of x, /x could be 


as large as 107! in some conditions, and as small as 10~* in other conditions. 


Recently, Qin et al. (2013) proposed that the shock acceleration strength makes important 
contributions to the reservoir phenomenon, particularly in low-energy SEPs. In their simulations, 
the reservoir phenomenon is reproduced under a variety of conditions of shock acceleration 


strength and perpendicular diffusion. In this paper, as a continuation of Qin et al. (2013), we study 
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the property of SEP spectra in the decay phase. We compute the time profiles of SEP flux which 
are accelerated by interplanetary shock. In section 2 we describe the SEP transport model and the 


shock model. In Section 3 we show the simulation results. In Section 4 we summary our results. 


2. MODEL 


In this work, we model the transport of SEPs following previous research(e.g., Qin et al. 
2006; Zhang et al. 2009; Dróge et al. 2010; He et al. 2011; Zuo et al. 2011; Wang et al. 2012; Qin 
et al. 2013; Zuo et al. 2013; Wang et al. 2014). A three-dimensional focused transport equation is 
written as (Skilling 1971; Schlickeiser 2002; Qin et al. 2006; Zhang et al. 2009) 


af ^ à (, ðf 
— =V. -V = weVv D 
2: (k,- Vf) (wb ev | feel u x 
Ho 

+p 1 H (v.v -bi : wv) +1285: vye A 
l-wy v AA mE 

_ ane V.ysw — : VV“ || — 1 

2 | pH vm -388: vv J 3: » 


where f(x, u, p.t) is the gyrophase-averaged distribution function; x is the position in a 
non-rotating heliographic coordinate system; f is the time; u, p, and v are the particle pitch-angle 
cosine, momentum, and speed, respectively, in the solar wind frame; b is a unit vector along the 
local magnetic field; V^" = V?" f is the solar wind velocity; and L is the magnetic focusing length 
given by L = (2 VinB) with Bo being the magnitude of the background magnetic field. The 
IMF is set as the Parker field model, and the solar wind speed is 400 km/s. This equation includes 
many important particle transport effects such as particle streaming along the field line, adiabatic 
cooling in the expanding solar wind, magnetic focusing in the diverging IMF, and the diffusion 


coefficients parallel and perpendicular to the IMF. 


The pitch angle diffusion coefficient model is set as (Beeck & Wibberenz 1986; Qin et al. 
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Dy = Doup?” lul! +h} (1-10), Q) 
where the constant Do controls the magnetic field fluctuation level. The constant q is chosen as 
5/3 for a Kolmogorov spectrum type of the power spectra density of magnetic field turbulence 


in the inertial range. Furthermore, h = 0.01 is chosen for the non-linear effect of pitch-angle 


diffusion at u = 0 in the solar wind (Qin & Shalchi 2009, 2014). 


The parallel mean free path (MFP) 4j can be written as (Jokipii 1966; Hasselmann 1968; Earl 


1974) 
3v (thd -,2y 


A = 5 
8 J-a Dy 


du, (3) 
and the parallel diffusion coefficient xj can be written as xy = v44/3. 


The relation of the particle momentum and the perpendicular diffusion coefficient is set as 


(Potgieter & Moraal 1985; Zhang 1999) 


«(rase (z) - 66) e 


where B, is magnetic field strength at the Earth, B is the magnetic field strength at the location of 
particle, p is particle momentum, and a is set to 1/3. Different perpendicular diffusion coefficients 
could be obtained by altering xo. Note that we use this ad-hoc model for the purpose of simplicity, 
the parameters, e.g., œ could be set as other values (Zhang 1999). However, the variation of these 
parameters would not qualitatively change the results in this paper. There are some more complete 
models that are developed to describe the particle diffusion in magnetic turbulence, such as the 
nonlinear guiding center theory (Matthaeus et al. 2003; Shalchi et al. 2004, 2010; Qin & Zhang 
2014). 


We use a time-backward Markov stochastic process method to solve the transport equation 
(1). The detail of method can be found in Zhang (1999) and Qin et al. (2006). The particle 


injection on the shock is specified by boundary values. The boundary condition is chosen as 
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following form (Kallenrode & Wibberenz 1997; Kallenrode 2001; Wang et al. 2012; Qin et al. 
2013) 


far, 0, 9, p. t) = a-ó(r — ust) S(r,6, 9, p): p ? - £(0,) 


0, 
S(r0,p,p) = C -exp |- E 
&0,9) = BONS? P 


0 otherwise, 


where the particle are injected at r = ust, v, is shock speed. ust = ro +n: Ar, with n = 0, 1,2 -+ - no. 
Ar is space interval between two ‘fresh’ injections, ry = 0.05 AU is inner boundary. r, is set to 
0.05. r is distance between sun and shock. ¢ is the angle between the center of shock and the 
point at the shock front where the particles injected. The shock acceleration strength is set as S for 
specifying the particles ejection. It changes with a power law in radial distance and exponential 
towards the flank of shock. £ determines the spatial scale of shock front. $, is the half width of the 


shock. 


3. RESULTS 


The parameters used are listed in table 1, unless otherwise stated. Note that the IMF is set 
as Parker spiral, and the disturbances of IMF behind the shock are ignored. The particle energy 
channels are chosen as 5 MeV, 10 MeV, 20 MeV, 40 MeV and 80 MeV. The parallel mean free 
path depends on the momentum A, ~ p!?. According to Qin et al. (2013), the x, /«j is set as 
0.1 in the ecliptic at 1 AU. Because the shock acceleration efficient decreases as the particle 
energy increasing, the acceleration strength parameters also change with the momentum: a ~ p9?, 


$e ~ p 99. The observers are located in the ecliptic at 1 AU. 
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3.1. Temporal Invariance in the Spectra 


In Figure 1 and 2, we plotted the fluxes of different energy channels in the cases with and 
without adiabatic cooling. In order to check the temporal properties of SEP spectra in the decay 
phase, we normalize different energy fluxes, so that the fluxes have similar values soon after all of 
them have reached peaks. The panels show the normalized fluxes observed in the ecliptic at 1 AU, 
but at different longitudes E60, E20, W20 and W60. The notations E60, E20, W20, W60 are short 
for East 60^, East 20°, West 20° and West 60°, respectively. East/West means the location of the 
Observer is east/west relative to the center of shock. The vertical lines indicate the shocks passage 


of the observers. 


In Figure 1, the adiabatic cooling effect is included in SEP propagation process. In the 
decay phase of SEP events, shock acceleration strength, adiabatic cooling, parallel diffusion, and 
perpendicular diffusion are the major factors to influence the flux behavior. In the four panels, the 
fluxes at all energies follow a similar trend, then the fluxes scatter slowly as time goes by. This is 
called temporal invariance in the spectra of gradual SEP event. In the E60 event, the fluxes at all 
energies start to follow a similar trend about one day before the shock passage of 1 AU. In other 
words, the onset time of the temporal invariance is earlier than the time of the shock passage of 
the observer. In the E20 and W20 events, however, the onset time of temporal invariance is close 
to the shock passage of the observers. In the W60 event, furthermore, the temporal invariance 


starts the latest, and actually it starts two days latter than the shock arrival. 


In Figure 2, the adiabatic cooling is not included in the SEP propagation process. Without 
adiabatic cooling, shock acceleration strength, parallel and perpendicular diffusion are the major 
factors in the decay phase. Due to the different diffusion coefficients and shock acceleration 
strength for different energy particles, the fluxes decay with different ratios consequently. With 
higher energies, the fluxes decay much faster. In these cases, the temporal invariance does not 


exist in the decay phase. 


199 


edis 


Comparing Figure 1 and Figure 2, the fluxes decrease much faster with adiabatic cooling. 
Because of adiabatic energy loss, particles have less energy when they are observed than that 
when they are released in the sources. In addition, since the source spectrum index is negative, 
the fluxes are lower with higher energies, so that the adiabatic cooling effect makes the SEP flux 
decreasing as time passes by. To sum up, the temporal invariance in the spectra results from the 


adiabatic cooling effect. 


3.2. Spatial Invariance in the Spectra 


In Figure 3, the SEP fluxes are shown for three observers located at different longitudes, E20, 
W20, and W60. The upper panel shows the 5 MeV proton fluxes observed by the observers. We 
set two typical time intervals, interval A from 1.3 days to 1.5 days in rising phase and interval 
B from 6.9 days to 7.1 days in decay phase. In order to study the spatial variance in different 
phases, in lower left and right panels of Figure 3, we plot the energy spectra observed in different 
longitudes in interval A and interval B, respectively. During the interval A, the spectra are different 
among the three observers. However, during the interval B, spectra are almost the same among the 
three observers. This phenomenon, which is named spatial invariance in the spectra by Reames 


et al. (1997), results from the reservoir effect in different energy channels. 


Because the shock is a moving source in the interplanetary space, the peak intensity of SEP 
flux is mainly determined by the shock acceleration strength and parallel MFP. In upper panel of 
Figure 3, at the peak time of flux for W20 (W60), the flux for W20 (W60) is close to that for E20. 
Furthermore, the SEP fluxes decay in a similar ratio because of the effect of adiabatic cooling. 
At the same time, the latitudinal gradient in the SEP fluxes is further reduced because of the 
effect of perpendicular diffusion. However, in other simulations with different shock acceleration 
strength and parallel MFP (not shown here), if at the peak time of flux for W20 (W60), the flux 


for W20 (W60) is significantly different than that for E20, the reservoir phenomenon can not form 
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in normal diffusion coefficients. As a result, shock acceleration strength, parallel diffusion, and 
adiabatic cooling are the main factors in forming the reservoir phenomenon, and perpendicular 


diffusion is a secondary factor. 


3.3. Invariant Spectra Region 


There are some important characteristics in the invariant spectra region from our simulations 
(Figure 1). If the observer is located at the eastern flank of the shock, the onset time of invariant 
spectra is earlier than the shock arrival. But if the observer is located near the central flank of 
shock, the spectra invariance begins approximately at the shock passages. Finally, if the observer 
is located at the western flank of shock, the onset time of invariant spectra are much later than the 


shock arrivals. From these results, we can better understand the invariant spectra region. 


Figure 4 shows the invariant spectra region. In the picture, the green line is plotted by 
Reames et al. (1997), and the red line is new in this work. According to Reames et al. (1997), the 
left side of the green line is the invariant spectra region, with the assumption that particles are 
quasi-trapped in the region behind the ICME, and the SEP fluxes gradually decrease as a result 
of parallel diffusion and adiabatic deceleration mechanisms, and in addition, there are also some 
leakage of energetic particles from ICME to eastern side of the upstream shock. In this sense, 
ICMEs play the pivotal role in the decay phase of fluxes. As a result, the invariant spectra region 
is determined by ICMEs’ propagation path plus some eastern side of upstream region. However, 
we suppose the invariant spectra region could be in the left side of the red line instead. In our 
simulations, ICME is not included, but in the propagation process perpendicular diffusion is 
included to reduce the spatial gradient in the fluxes, and adiabatic cooling is included to reduce 
the temporal variance. As the simulation results showed above, the spectra spatial and temporal 
invariance could result from the effects of shock acceleration strength, adiabatic cooling, and 


perpendicular diffusion. In this sense, it is possible that the invariant spectra region is not confined 
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by the ICMEs’ propagation path. Instead, the invariant spectra region could be confined by the 
interplanetary shock, but the region expands faster (slower) than the shock at the eastern (western) 


flank, respectively. 


4. DISCUSSION AND CONCLUSIONS 


We have studied interplanetary shock accelerated SEPs propagation in three-dimensional 
IMF. The spectra observed by different observers are calculated, and the spatial and temporal 


invariance in the spectra are reproduced in the simulations. The following are our major findings. 


The adiabatic cooling effect is the key factor for forming temporal invariance in the spectra. 
By including the adiabatic cooling, for different energy channels, the flux decay ratios are almost 
the same. The temporal invariance results from the fact that all energy particles decay as the same 
ratios because of adiabatic cooling effect. At the eastern flank of the shock, the onset time of the 
spectra invariance is earlier than the shock arrival. For the central cases, however, the onset time 
of the spectra invariance is close to the time of shock arrival. At the western flank of the shock, 
finally, the onset time of the spectra invariance is later than shock arrival. In addition, the fluxes 
decay much faster in the cases with adiabatic cooling. Without adiabatic cooling, the decay phase 
of SEP fluxes are dominated by shock acceleration strength, parallel diffusion, and perpendicular 
diffusion, which are all varying with particles’ energies. Therefore, the temporal invariance does 


not exist without adiabatic cooling. 


Shock acceleration strength, parallel diffusion, adiabatic cooling, and perpendicular diffusion 
are four important factors in forming the spatial invariance, which is in reservoir phenomenon in 
different energy channels. Shock acceleration strength parameters œ and ¢, are set to 2 and 15° 
for 5 MeV protons in our simulations, respectively. And these parameters also change with the 


momentum: a ~ p9?, de ~ p 93, because the shock acceleration strength decreases with higher 
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energy particles. Among the four factors, shock acceleration strength, parallel diffusion, and 
adiabatic cooling are the main factors in forming the reservoir phenomenon, and perpendicular 
diffusion is a secondary one. This conclusion is derived based on our simulations, and it is also 
consistent with the observations. In Reames et al. (1997) and Reames (2013), a gradual SEP event 
was detected by Helios 1, Helios 2, and IMP 8 during March 1, 1979 to March 11, 1979. In 
this event, the reservoir phenomenon appeared, and the onset time of SEPs observed by different 
spacecraft are very close because of the effect of perpendicular diffusion. At the peak time of flux 
observed by Helios 2 (IMP 8), the flux observed by Helios 2 (IMP 8) is close to that observed 
by Helios 1. The importance of the peak intensity of SEP flux observed by Helios 2 and IMP 8 
in forming the reservoir phenomenon is also noticed by Reames (2013), however, the reservoir 
phenomenon is explained as a result of the disturbance of IMF caused ICMEs. In our model, 
the peak of flux is mainly determined by shock acceleration strength and parallel diffusion. 
Furthermore, the SEP fluxes decay as a similar ratio because of the effect of adiabatic cooling. At 
the same time, the latitudinal gradient in the SEP fluxes is further reduced because of the effect of 
perpendicular diffusion. Finally, according to our model, the reservoir phenomenon appeared in 
this SEP event with the effects of shock acceleration strength, parallel diffusion, adiabatic cooling, 
and perpendicular diffusion. Observationally, shock acceleration strength, diffusion coefficients, 
and adiabatic cooling change significantly in different SEP events (Kallenrode 1996, 1997). As a 
result, the reservoir phenomenon can only form in some gradual SEP events with those controlling 


effect parameters in appropriate values. 


Based on our simulations, a new invariant region is proposed. The new region is different 
from the one proposed by Reames et al. (1997). There are two important characteristics in our new 
region. First, if the observer is located at the eastern (western) flank of the shock, the onset time of 
temporal invariant in the spectra is earlier (later) than the shock arrival, respectively. Second, the 
spatial invariance in the spectra can also be formed without ICMEs. These two characteristics are 


supported by observations, but are difficult to be explained in the previous model Reames et al. 
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(1997). 


In our model, we ignore the disturbance of the IMF caused by ICME for the simplicity. In 
principle, the disturbance in the magnetic field can help particles redistribute in space. In future 
work, we intend to include a realistic three-dimensional ICME shock, so that the SEP acceleration 


and transport in the heliosphere can be investigated more precisely. 
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Fig. 1.— Comparison of different energy protons observed by observers in the 1 AU ecliptic at 
different longitudes. The fluxes are normalized so their values are similar after all of them reach 
peaks. The vertical lines indicate the shocks passage of 1 AU. The adiabatic cooling is included in 


simulations. 
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Fig. 2.— Same as 1 except that the adiabatic cooling is not included. 
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Fig. 3.— In the upper panel, comparison of 5 MeV protons flux observed by the observers in 1 


AU ecliptic at different longitudes. The lower left and right panels show the spectra observed at 


different longitudes during time interval A and B, respectively. The vertical dashed lines indicate 


the shock passage of 1 AU. 
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Fig. 4.— The green line indicates the original spectra invariant region proposed by Reames et al. 


(1997). The red line indicates a new spectra invariant region based on the simulation results in this 


paper. 
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Table 1: Model Parameters Used in the Calculations. 


Parameter | Physical meaning Value 
yer solar wind speed 400 km/s 
Vs shock speed 870 km/s 
ds shock width 60° 

a shock strength parameter 25 

be shock strength parameter 15°? 

y injection spectrum 5.5 

B, magnetic field strength at the Earth | 5 nT 

Aj particle radial mean free path 0.2 AU* 
Kı perpendicular diffusion coeffient 0.1 x «qd 
po inner boundary 0.05 AU 
Fpi outer boundary 50 AU 


“for 5 MeV protons. 
>for 5 MeV protons. 
“for 5 MeV protons in the ecliptic at 1 AU. 


“for 5 MeV protons in the ecliptic at 1 AU. 


